Linking the Supermassive Black Hole Growth with 
the Megamaser Emission 
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Abstract. High-resolution observations of the central few 100 pc of the galactic nuclear 
environments remain prohibitive for large statistical samples, which are crucial for tracing the 
, links between central black hole formation, galaxy formation and AGN activity over cosmic 

time. With this contribution, we present novel ways of connecting the physics of black hole 
accretion with its immediate environs via a new quantitative evaluation of the degree to which 
pG ■ the strength and spatial configuration of the water maser emission is associated with the nuclear 

Oh' nebular galactic activity. We discuss possible evolutionary/causal connections between these 

two types of emission, together with criteria that could dramatically enhance our search for 
5^ ■ mega-maser systems in nearby galactic centers. These are timely results given the interest in 

' combining high-resolution observations with extremely large optical telescopes and large arrays 

C^j , that start to conquer the sub-millimeter window. 
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£S| , 1. Introduction 

Mega-masers are remarkable natural phenomena. Mapping of the H2O maser emission with 
Ln 2 o > IOLq in a disk configuration has provided the first direct evidence in an active 
galactic nucleus (AGN) for the existence of a thin Keplerian accretion disk with turbulence 
|15j . highly compelling evidence for the existence of a massive black hole (BH) [281 I27| . as 
well as a cosmic distance determination of extremely high precision \19\ 120] . To date, such 
observations and associated calculations have been performed on only less than a handful such 
systems |31|, [27] . Identifying and studying a large sample of the mega-masers in this particularly 
rare disk configuration is of crucial relevance for our progress in understanding the accretion 
processes in galaxy centers, and more importantly, for constraining the proposed cosmological 
models and the nature of Dark Energy. Unfortunately, good candidates are extremely few; the 
success rate in maser detection remains a mere ~ 3 — 5%, out of which only ~ 40% appear to 
originate in disks H QH HU E] . 

To efficiently search for new powerful mega-maser disks, we need a good understanding of the 
special physical characteristics that facilitate their production in galaxy centers. There is some 
evidence that mega-masers may be associated with the molecular disk or torus that surrounds 
and (partially) obscures an actively accreting massive black hole (Mbh ~ 10 7 M sun ) harbored 
by a galactic nucleus, i.e., an AGN [281 IS]- However, mega-masers are also detected in face-on, 
unobscured disks or elliptical galaxies where no molecular torus is expected, and even along jets 
protruding long distances from the nucleus [12l|T3]. So, either masers are not always produced 
in association with an AGN torus, or these particular cases are unrecognized AGN. Indeed, for 
the majority of the nearby active galaxies, an AGN classification remains ambiguous; in these 



sources, an actively accreting black hole might be either obscured by large amounts of dust in the 
surrounding star- forming regions, or simply hidden in a mix of other ionization sources (shocks, 
turbulence, etc.) or host galaxy light [21] . Thus, if mega-maser production is genuinely related 
to black hole accretion, the true detection rate of disk mega-maser emission remains unknown 
as long as the AGN census is incomplete. To shed light on these matters, we need to investigate 
all possible connections between the H2O emission and their host properties, for a wide variety 
of galactic nuclear activity. We only now have the number statistics necessary to perform such 
a study. 

We introduce here an in-depth multi-parameter characterization of the optical properties of 
the largest available sample of 146 galaxy centers hosting H2O masers together with a ~ 20 times 
larger control sample of galaxies in which maser detection failed. This effort complements the 
very recent work by Zhu et al. [3B], where they investigate the optical properties of the maser 
host galaxies and find that the maser detection rate increases for higher optical luminosity, larger 
velocity dispersion, higher extinction and higher [O III] A5007 line luminosity. Our analysis adds 
a novel comparative investigation of the optical host and nuclear features separated per spectral 
types that correlate with the strength and type of nuclear galactic activity. If, as previously 
suggested [91 [101 [52~| 133] . there is an evolutionary sequence in which galaxies transform from 
star-forming via AGN to quiescence, accretion onto the central BH should happen as early as 
the star formation (or H II) phase, and should be active in Transition Objects (or Ts). Thus, 
if maser activity is related to accretion, we should look for (mega)masers in H n's and Ts as 
well. Restricting our searches to Seyferts (Ss) and Low Ionization Nuclear Emission Regions (or 
LINERs; Ls), as it appears to be the case with current efforts, could overlook a good portion, 
or even the lion's share of mega-masers. 

2. The Sample Galaxies &; Parameters Involved 

The Mega-maser Cosmology Project (MCP; |31] IB]) makes publicly available the entire sample 
of galaxies detected in H2O by May 2011. The data is presented in the form of an atlas that 
comprises sky positions, recession velocities, the maser spectra, and the corresponding discovery 
reference, for 146 such systems. The MCP atlas also includes the list of all > 3300 galaxies that 
have been surveyed for H2O maser emission, along with their Green Bank Telescope (GBT) 
spectra, recession velocities, the sensitivity of each observation, and the corresponding source 
brightness temperature. These data provide for the first time a sufficiently large sample of 
maser galaxies and a control sample of non-detections whose properties can be compared at a 
statistically significant level. 

Optical spectral and photometric measurements for these objects are obtained mainly from 
wide sky surveys; for the purpose of the analysis presented here, we exploited the nearby galaxy 
sample of the seventh data release of the SDSS [I]. Measurements of the host and nuclear nebular 
emission of the SDSS galaxies are drawn from the catalog built by the MPA/JHU collaboration 
[7]. For these objects, the emission-line component is measured after it is separated from the 
host stellar emission, and subtracted from the total galaxy spectrum based on fits of stellar 
population synthesis templates. Black Hole masses (Mbh) and corresponding accretion rates 
dM/dt (measured via the Eddington ratio L/Lsdd) are obtained based on estimates of the 
stellar velocity dispersion (c*), that are available via careful spectral fitting of the host stellar 
light component in galaxy spectra. To relate the maser and the central BH accretion activity 
to the host properties, we employ stellar masses and the -D n 4000 break index or H<5yi Balmer 
absorption-line index as proxies for the age of the associated stellar population, as calculated 
and presented by Kauffman et al. |22| [23] . 

Cross-matches of the maser and control galaxy lists with these SDSS-based catalogs uncover 
79 maser and ~ 1600 control galaxies with photometric information, among which high quality 
spectroscopic measurements (> 2a confidence) are available for 46 and 1181 maser and non- 



maser systems respectively. We present here only a comparison between the spectroscopic 
characteristics. It may be important to note that considering only the spectroscopic samples 
the 3% maser detection rate remains in place. 

3. A statistical comparison of the properties of maser and non-maser galaxies 

Optical emission line ratios have often proven to be useful in identifying the dominant ionization 
mechanism or in delineating differences between various types of galaxy nuclei [21 [35j I25j . A 
comparison of the emission-line activity in maser and non-maser galaxies is presented in Figures 
[T]and [2] where we also show how the distributions of the line ratios compare for the two samples. 
To classify the galaxies into Ss, Ls, Ts, and H us, we adopt the criteria of Kewley et al. [25J. 

With the exception of the [O III/H/3] ratio, which is clearly larger for the maser sample, 
there is no obvious statistical difference between the maser and non-maser (control) galaxies. 
The difference in [O III/H/3] could simply be a consequence of the fact that the bulk of the masers 
have been detected in surveys of Seyferts, whose definition requires high such ratios; however, 
when only Seyfert maser and control galaxies are compared, the difference is still statistically 
significant (KS pro b = 0.053), implying that it reflects an intrinsic relation to maser activity. 

The maser detection rate among Seyferts is high (~ 8%), however, it is unclear whether this 
is because the disk- maser emission is genuinely connected to this type of nuclear activity (i.e., 
the BH accretion process) or is simply a consequence of the same survey bias. The maser disks 
are usually mega-masers, but not all mega-masers are disks. The disks have only been found 
among the Seyferts, however, mega-masers appear to cover the whole spectral spectrum. The 
maser detection rate is clearly non-zero for non-AGN types, in fact, there is an equal probability 
to find masers among Ts as in Ls (~1.6%), and twice as likely to find them in H us (~ 3%). 
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Figure 1. Line-diagnostic diagram for all maser and control galaxies, 
with a comparison of histograms of the involved line ratios. The non- 
maser galaxies are plotted with grey symbols, the masers in black, the 
mega-masers (Lh 2 o > 10L S un) are shown as magenta squares while the 
disks among them are the blue circles. The lines in the bottom panels 
are semi-empirical fits to the distributions of objects in these diagrams. 
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Figure 2. Comparison of the 
distributions of [O IIfJ/H/3 line ra- 
tios for maser and control galaxies. 
[O III]/H/3 is clearly higher in maser 
hosts. 



3.1. The host characteristics that differentiate between maser and non-maser galaxies 
Our comparisons of overall distributions of a variety of optical host properties of maser and 
non-maser galaxies support previous findings that: 1) maser galaxies show statistically higher 
levels of reddening and extinction, as measured by the Hce/H/3 Balmer decrements; only the 



maser galaxies show values of Ha/H/3 > 10; 2) maser galaxies are more massive (higher stellar 
mass M*), more luminous (in M r , not shown here), and harbor more massive BHs in their 
centers (as shown by the comparison in a*). We also show here for the first time that there is a 
possible connection between the maser activity and the age of the stellar population associated 
with their host galactic nuclei; the maser galaxies exhibit younger stellar populations (weaker 
-D n 4000, suggesting an order of magnitude difference in age, in average). Figure [3] presents 
these results. 




Figure 3. Distributions of redshifts z, Balmer Decrements Ha/Hf3, stellar velocity dispersions a* , host 
stellar mass M* , and £>„4000 break index for maser galaxies (black lines) and the control sample (dashed 
gray lines). 

Given the apparent bias of the current maser surveys toward Seyfert galaxies, interpreting 
these global trends might benefit from looking into how these properties compare when maser 
and non-maser galaxies are separated per spectral types, as defined based on the line-diagnostic 
diagrams (Figure [1]). Figure H] illustrates the individual measurements of redshift, Balmer 
decrements, stellar masses, and L>„4000 parameters, for maser and non-maser galaxies, separated 
by their respective spectral types H II, T, S, and L. This comparison reveals some features worth 
mentioning: i) H n's with masers are only at very low z, and only if they are dwarfish, i.e., with 
very low stellar masses; ii) Ts associated with maser emission are more nearby, less massive and 
clearly younger in their stellar populations than their counterparts without maser activity; iii) 
interestingly, Ss with and without masers don't separate much from each other in these host 
properties; iv) it is only the high-z Ls with the youngest stellar populations that host masers, 
which are all mega- masers. Interestingly, the mega-maser disks are not among the galaxies with 
the youngest stellar populations, or among those with extremely large Balmer decrements. It 
is probably surprising that the only H n hosting mega-maser emission has the highest Balmer 
decrement in the whole maser galaxy sample; nevertheless, this is not a mega-maser disk. 

3.2. The nuclear features that differentiate between maser and non-maser galaxies 
A very similar type of comparison for measures that pertain more directly to the process of 
BH growth in these systems is presented in Figure [5j Here we show individual measurements 
in 1) stellar velocity dispersion a*, as a measure of the BH mass; 2) [O III] line luminosity, 
L[OIII], which is corrected for reddening and extinction based on a r oc A -07 attenuation law 
[8] and the corresponding Balmer decrements; 3) a measure of the Eddington ratio L/L e ^, 
where L = L\> \ = 600 x L[OIII] was used for the bolometric correction [181 121] and L e dd is 
the Eddington luminosity corresponding to the Mbh values obtained via Mbh — o~* relation of 
Graham et al. [T3]; and 4) the [S II] A 6716/A6731 line ratio as a proxy for the density n e of 




Figure 4. Host properties compared as a 
function of the different optical spectral types H n, 
T, S, and L. Symbol types are as in Figure [1] 
Large open squares indicate averages per spectral 
type. 



Figure 5. Nuclear properties compared as a 
function of the different optical spectral types H II, 
T, S, and L. Symbol types are as in Figure [TJ 
Large open squares indicate averages per spectral 
type. 



the line-emitting gas [30J, which may also be viewed as a gauge of the availability of fuel for 
accretion and/or masing activity. 

We find here that the previously claimed trends for higher BH masses and L[OIII] in maser 
galaxies [3B] are correct not only for the Seyfert hosts, but for galaxies of all spectral types. 
While it seems especially true that among Ls, maser activity is detected only in systems with 
the highest L[OIII] and the most massive BHs, Seyferts with and without masers exhibit a 
relatively weak statistical trend in this direction. Note also that galaxies with very massive BHs 
(a* > 200 km/s) do not host masers; physical conditions in the vicinity of such systems may 
not be stable enough to sustain strong mega-maser emission. 

Additionally, our comparison, reveals some new and potentially very useful constraints that 
could be posed to future maser surveys, particularly when searching for mega-maser disks: 

1) the combination of Mbh and L[OIII] used to calculate L/L e ^ exposes a potentially narrow 
range in the accretion rate of nuclei associated with maser emission, which in average is higher 
than in the control sample, for all the different spectral types compared here; it seems that the 



maser activity, and especially the mega-maser disk-like emission, is associated to accretion onto 
BHs characterized by logLjL^A ~ —2. Incidnetally, this value also appears to correspond 
to a change in the accretion flow geometry or mode, from efficient (i.e., Shakura-Sunyaev 
standard accretion disk/corona, or quasar- like, [M]) to inefficient (i.e., the radiatively inefficient, 
or advection dominated accretion flow; RIAF/ADAF; e.g., |29| ) with decreasing L/L e ^d [TOlfTT]. 
It may thus be possible that the disk mega-maser emission is related to a certain (short) phase 
in the galaxy-AGN evolution, reflected by this particular change. 

2) the comparison of the [S II] A 6716/A6731 line ratio unfolds a potential "goldilocks range" 
in n e for the maser galaxies, spanning only a range of ~ 100 — 300 cm -3 (corresponding to a range 
in [S II] A 6716/A6731 = 1.4 — 1). The galaxies with the highest and lowest n e in the surveyed 
sample do not host maser activity at the levels that can be currently detected. Because these 
measurements of n e reflect physical conditions of regions located way outside the masing areas, 
this parameter relates to the maser detection rate only indirectly; nevertheless, this connection 
is not less essential than that between L[OIII] and the maser activity, which suffers from similar 
caveats. 

4. Summary 

We presented here some results of the first multi-parameter comparative analysis of the largest 
sample available of SDSS-based galaxies with and without maser activity in their centers. Besides 
confirming previously claimed tendency for high L[OIII] in maser galaxies, we found new criteria 
for more efficient targeting of maser galaxies, and more importantly, of hosts of mega-maser 
disks, based on their optical spectral properties. Our study suggests that mega-masers could be 
found in emission-line galaxies that are not necessarily classified spectroscopically as Seyferts, 
and that more appropriate survey designs should consider, at least to start with, a cross-match 
of narrow ranges in their L/L^ values and the [S II] A 6716/A6731 line ratios, as described 
above. Considering a survey of all galaxies included in this study with —2 < logL / L c ^ < and 
[S II] A 6716/A6731 = 1 — 1.4, the maser detection rate is boosted by a factor of 4, to ~ 12%, 
while the maser disk detection rate increases to > 25%. More details on this analysis, along 
with a comprehensive Principal Analysis that encompasses a wider variety of parameters than 
those included in this discussion will be available in Constantin et al [TT] . 
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